REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

23-05-2006  Journal  Article 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

Electronic  Structure  Studies  of  Tetrazolium-Based  Ionic  Liquids 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Deborah  D.  Zorn  &  Mark  S.  Gordon  (Iowa  State  Univ.);  Jerry  A.  Boatz  (AFRL/PRSP) 

5d.  PROJECT  NUMBER 

23030423 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRSP 

10  E.  Saturn  Boulevard 

Edwards  AFB  CA  93524-7680 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-PR-ED-JA-2006-289 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

Air  Force  Research  Laboratory  (AFMC) 
AFRL/PRS 
5  Pollux  Drive 


11.  SPONSOR/MONITOR’S 
NUMBER(S) 


Edwards  AFB  CA  93524-7048 


AFRL-PR-ED-JA-2006-289 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited  (AFRL-ERS-PAS-2006-027) 


13.  SUPPLEMENTARY  NOTES 

Published  in  Journal  of  Physical  Chemistry  B  2006,  110,  11110-11119 


14.  ABSTRACT 

New  energetic  ionic  liquids  are  investigated  as  potential  high  energy  density  materials.  Ionic  liquids  are  composed  of  large,  charge- 
diffuse  cations,  coupled  with  various  (usually  oxygen  containing)  anions.  In  this  work,  calculations  have  been  performed  on  the  tetrazolium 
cation  with  a  variety  of  substituents.  Density  functional  theory  (DFT)  with  the  B3LYP  functional,  using  the  6-3 1  lG(d,p)  basis  set  was  used  to 
optimize  geometries.  Improved  treatment  of  dynamic  electron  correlation  was  obtained  using  second  order  perturbation  theory  (MP2).  Heats  of 
formation  of  the  cation  with  different  substituent  groups  were  calculated  using  isodesmic  reactions  and  Gaussian-2  calculations  on  the  reactants. 
The  cation  was  paired  with  oxygen  rich  anions  CI04',  N03",  or  N(N02)2’  and  those  structures  were  optimized  using  both  DFT  and  MP2.  The 
reaction  pathway  for  proton  transfer  from  the  cation  to  the  anion  was  investigated. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 
PERSON 

Dr.  Scott  A.  Shackelford 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

A 

11 

19b.  TELEPHONE  NUMBER 

(include  area  code) 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


11110 


J.  Phys.  Chem.  B  2006,  110 ,  11110-11119 


Electronic  Structure  Studies  of  Tetrazolium-Based  Ionic  Liquids 

Deborah  D.  Zorn,1^  Jerry  A.  Boatz,* *  and  Mark  S.  Gordon*^ 

Department  of  Chemistry,  Iowa  State  University,  Ames,  Iowa  50011,  and  Air  Force  Research  Laboratory, 

Space  and  Missile  Propulsion  Division,  AFRL/PRS,  10  East  Saturn  Boulevard,  Edwards  AFB,  California  93524 

Received:  February  9,  2006;  In  Final  Form:  April  3,  2006 


New  energetic  ionic  liquids  are  investigated  as  potential  high  energy  density  materials.  Ionic  liquids  are 
composed  of  large,  charge-diffuse  cations,  coupled  with  various  (usually  oxygen  containing)  anions.  In  this 
work,  calculations  have  been  performed  on  the  tetrazolium  cation  with  a  variety  of  substituents.  Density 
functional  theory  (DFT)  with  the  B3LYP  functional,  using  the  6-311G(d,p)  basis  set  was  used  to  optimize 
geometries.  Improved  treatment  of  dynamic  electron  correlation  was  obtained  using  second-order  perturbation 
theory  (MP2).  Heats  of  formation  of  the  cation  with  different  substituent  groups  were  calculated  using  isodesmic 
reactions  and  Gaus sian-2  calculations  on  the  reactants.  The  cation  was  paired  with  oxygen  rich  anions  ClOzT, 
N03-  or  N(NC>2)2~  and  those  structures  were  optimized  using  both  DFT  and  MP2.  The  reaction  pathway  for 
proton  transfer  from  the  cation  to  the  anion  was  investigated. 


Introduction 

There  is  considerable  current  interest  in  ionic  liquids  as 
solvents.  No  other  class  of  solvents  offers  the  versatility  that 
ionic  liquids  do:  They  are  typically  thermally  stable  and  have 
negligibly  low  vapor  pressure,  high  density,  and  large  liquid 
ranges  up  to  400  °C.  Because  of  these  unique  properties,  ionic 
liquids  may  be  used  as  electrolytes  for  batteries,  extraction 
media,  and  catalyst  carriers.  The  popularity  of  ionic  liquids  has 
also  been  spurred  on  by  their  classification  as  “green”,  due  to 
their  negligible  vapor  pressure,  thereby  decreasing  levels  of 
volatile  organic  carbons  in  the  environment.1 

l-Butyl-3-methyl-imidazolium  (BMIM)  is  one  of  the  most 
widely  used  and  best  understood  ionic  liquid  cations,  so  when 
one  is  beginning  to  design  energetic  ionic  liquids,  this  cation  is 
a  natural  place  to  start.  The  unsubstituted  imidazolium  cation 
is  shown  in  Figure  la.  The  energy  content  of  the  ionic  liquid 
cation  needs  to  be  raised  for  high  energy  applications,  and  this 
can  be  done  in  several  ways. 

First,  replacing  the  imidazolium  ring  with  a  more  nitrogen 
rich  ring  can  increase  the  energy  content  of  the  ionic  liquid. 
Second,  the  hydrogen  or  alkyl  side  chains  can  be  replaced  with 
high  energy  groups,  such  as  — CN,  —  NH2,  —  N3,  and  —  NO2. 
Finally,  the  anion  chosen  should  be  oxygen  rich  to  serve  as  an 
oxidizer.  Interesting  anions  include  nitrate,  perchlorate,  and 
dinitramide  anions  (Figure  2).  Two  promising  nitrogen-rich 
candidates  are  the  triazolium  (Figure  lb)  and  the  tetrazolium 
(Figure  lc)  cations  with  three  and  four  nitrogens  respectively, 
in  the  ring.  These  two  cations  could  present  a  useful  balance 
between  exothermicity  and  thermal  stability.  A  thorough  study 
of  the  1,2,4-triazolium  cation  has  recently  been  published.2  The 
focus  of  the  present  study  is  on  the  electronic  structure  of  the 
tetrazolium  cation. 

The  energetic  tetrazolium  cations  that  have  been  synthesized 
include  1 -amino-4, 5-dimethyltetrazolium, 3  2-amino-4,5-dimeth- 
yltetrazolium,3  2,4,5-trimethyltetrazolium,4  and  l,5-diamino-4//- 


*  Corresponding  author. 

t  Department  of  Chemistry,  Iowa  State  University. 

*  Air  Force  Research  Laboratory,  Propulsion  Sciences  and  Advanced 

Concepts  Division,  AFRL/PRS. 


a.) 


H — n;  n — H 


b.) 

H 


H — XC— H 

W  J 

N — N 

\ 

H 


c.) 

H 

I 

/N1x 

H — C5  N2 

W  J! 

N4-N3 

\ 

H 


Figure  1.  Imidazolium  (a),  triazolium  (b),  and  tetrazolium  (c)  cations. 
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Figure  2.  Nitrate,  perchlorate  and  dinitramide  anions  (oxygen  =  red, 
nitrogen  =  blue,  and  chlorine  =  green). 

tetrazolium.5  Melting  points  for  these  cations  combined  with 
iodide,  nitrate,  and  perchlorate  anions  range  from  —59  to  +156 
°C  and  their  liquid  range  can  be  up  to  229  °C.  A  summary  of 
energetic  tetrazolium  cations  that  have  recently  appeared  in  the 
literature  and  their  melting  points  are  given  in  Table  l.3-5 

Synthesis  of  triazolium  cations  is  easier  than  that  of  tetrazo¬ 
lium  cations,  hence,  there  have  been  many  more  studies  of 
energetic  triazolium  than  tetrazolium  compounds.  Triazolium 
cations  have  been  successfully  substituted  with  azido,4  nitro,4 
and  amino  groups.3-6  The  successful  synthesis  of  the  triazolium 
cation  is  not  so  surprising  due  to  its  similarity  to  the  popular 
imidazolium  cation,  but  what  is  surprising  is  the  use  of  the 
triazole  ring  as  an  ionic  liquid  anion.  Ionic  liquids  typically 
contain  a  large  asymmetric  organic  cation,  which  causes  the 
ions  to  be  poorly  coordinated,  but  as  shown  by  Katritzky  et 
al.,7  the  anion  can  be  large  as  well.  An  ionic  liquid,  l-butyl-3- 
methylimidazolium  3,5-dinitro-l,2,4-triazolate,  containing  a  high 
energy  planar  anion,  has  been  synthesized  and  has  a  low  melting 
point  of  35  °C.  Delocalization  of  charge  on  the  anion  ring  is 
caused  by  the  nitro  substituent  groups,  which  also  serve  to 
increase  the  energy  of  the  compound. 

Interactions  in  ionic  liquids  are  more  complicated  than  those 
in  simpler  liquids,  making  them  more  difficult  to  understand 
on  a  molecular  level.  Theory  can  provide  an  excellent  tool  for 
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Figure  3.  Resonance  structures  of  isomers  I— IV. 

TABLE  1:  Energetic  Tetrazolium  Cations  and  Their  Melting  Points  (Tm),  Thermal  Decomposition  Temperatures  (Td),  and 


Glass  Transition  Temperature  (Tg) 

1  -amino-4, 5-dimethyltetrazolium  iodide" 

1  -amino-4, 5-dimethyltetrazolium  nitrate" 

1  -amino-4,5 -dimethyltetrazolium  perchlorate" 
2-amino-4, 5-dimethyltetrazolium  iodide" 
2-amino-4, 5-dimethyltetrazolium  nitrate" 
2-amino-4,5 -dimethyltetrazolium  perchlorate" 

2,4,5-trimethyltetrazolium  iodide^7 

2,4,5-trimethyltetrazolium  nitrate^7 

2.4.5- trimethyltetrazolium  perchlorate^7 

1 .5- diamino-4//— tetrazolium  perchlorate" 


Tm  121  °C 

Tg  -59  °C 

Td  170  °C 

rm  51  °c 

Td  182  °C 

Tm  124  °C 

Tm  94  °C 

Td  173  °C 

Tm  140  °C 

Td  238  °C 

Tm  156  °C 

Tm  94  °C 

Td  193  °C 

Tm  133  °C 

rd315°C 

Tm  125-130  °C 

"  From  Xue  et  al.3  b  From  Xue  et  al.4  c  From  Drake  et  al.5 


understanding  the  structure  and  dynamics  of  ionic  liquids.  There 
have  been  several  molecular  dynamics  (MD)  simulations  of 
ionic  liquids  based  on  imidazolium,  pyridinium  and  ammonium 
cations.  These  studies  modeled  bulk  properties,  such  as  melting 
points,8  diffusion,8-12  and  viscosity.8,9  Radial  distribution 
functions8,10,12-16  and  densities  have  also  been  calculated.8’10,15’17 
Dynamics  simulations  have  revealed  that  1 -ethyl-3 -methylim- 
idazolium  nitrate  has  diffusion  properties  similar  to  those  of  a 
supercooled  liquid.16  Systematic  MD  studies  of  l-alkyl-3meth- 
ylimidazolium  ionic  liquids  with  several  anions  have  been 
performed  using  various  classical  force  fields.11,17  Three  first 
principles  based  MD  studies  on  1,3-dimethylimidazolium 
chloride  have  been  performed  and  compared  to  classical  sim¬ 
ulations  and  experimental  neutron  scattering  experiments.12-14 
There  are  differences  in  the  local  structures  predicted  by  classical 
vs  first  principles  based  MD.  Efforts  to  improve  both  classical 
and  first  principles  based  MD  simulations  of  ionic  liquids  are 
ongoing  in  several  groups. 

The  solvent  properties  of  ionic  liquids  and  their  propensity 
to  be  synthesized  depend  greatly  on  their  acid— base  properties, 
so  the  ability  to  predict  these  properties  would  be  of  great  use. 
The  more  acidic  the  cation,  the  more  difficult  it  will  be  to 
protonate  the  neutral  ring  to  form  an  ionic  liquid.  Of  the  three 
neutral  rings,  imidazole,  triazole,  and  tetrazole,  the  weakest  base 
is  tetrazole.  Tetrazolium  is  therefore,  the  most  difficult  to 
synthesize.  Acidity  constants  of  triazole  and  tetrazole  have  been 
calculated  using  several  semiempirical  methods.19 

To  have  low  melting  ionic  liquids,  the  charge  of  at  least  one 
of  the  ions  must  be  delocalized.  Multiconfiguration  self- 
consistent  field  (MCSCF)  analysis  of  the  triazolium  cation  shows 
that  the  electrons  on  the  cation  are  shared  between  two  resonance 
structures.2  This  study  also  investigated  the  effects  of  more 
energetic  and  less  energetic  substituent  groups.  For  example, 
the  nitrile  group  was  proposed  as  a  better  substituent  for  high 
energy  applications  than  an  azide  group.  The  structures  of 
triazolium  dinitramide  systems  were  investigated  using  dimer 
pairs.  A  wide  variety  of  geometries  were  found  and  the  presence 


of  small  barriers  for  proton  transfer  from  the  cation  to  the  anion 
show  that  deprotonation  may  be  an  important  mechanism  in 
decomposition  of  triazolium-based  ionic  liquids. 

Computational  Methods 

Initial  structures  were  obtained  by  performing  gas-phase 
density  functional  theory  (DFT)  calculations  on  the  isolated  ions 
using  the  Becke  three  parameter  Lee— Yang— Parr  hybrid  func¬ 
tional  (B3LYP). 20,21  The  basis  sets  used  were  6-3  lG(d,p), 22-24 
and  6-311G(d,p).25  The  DFT  geometries  and  energies  are 
compared  to  second-order  Moller— Plesset  perturbation  theory 
(MP2)  calculations.26,27  A  MCSCF  population  analysis28  using 
Edmiston— Ruedenberg  type  localized  molecular  orbitals  (LMO)29 
was  carried  out  to  investigate  the  amount  of  electron  delocal¬ 
ization  in  the  ring.  The  orbitals  included  in  the  MCSCF  active 
space  are  both  Jt  orbitals,  their  corresponding  antibonds,  and 
the  lone  pair. 

Calculations  on  ion  pairs  provide  information  about  the 
fundamental  interactions  between  the  cation  and  the  anion.  The 
gas-phase  ion  interactions  can  provide  insight  into  the  bulk  liquid 
structure.  Dimer  pairs  were  optimized  using  DFT  and  MP2 
methods  and  the  6-3 1+G(d)23, 24,30,31  basis  set.  Hessians  (matrices 
of  energy  second  derivatives)  are  used  to  determine  whether 
stationary  points  are  minima  or  transition  states.  At  the  final 
MP2/6-311G(d,p)  geometries,  improved  relative  energies  were 
obtained  for  some  of  the  tetrazolium  cation  isomers,  using 
singles  and  doubles  coupled  cluster  theory  with  perturbative 
triples  (CCSD(T))32,33  with  the  6-31  lG(d,p)  and  cc-pVTZ34  basis 
sets. 

Intrinsic  reaction  coordinate  (IRC)  calculations  (also  referred 
to  as  the  minimum  energy  path)  were  used  to  connect  transition 
states  with  reactants  and  products.35,36  The  step  size  used  for 
the  IRC  calculations  was  0.05  (amu)1/2  bohr.  All  calculations 
were  done  with  GAMESS, 37,38  and  all  molecules  were  visualized 
with  MacMolPlt.39 
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Figure  4.  Open  chain  isomer,  azidoformidinium,  V. 


Figure  5.  MP2/6-31  lG(d,p)  minimum  energy  path  for  proton  transfer 
between  isomer  I  and  isomers  II,  III,  and  IV,  kcal/mol.  Structures  for 
the  three  transition  states  (I  — *  II,  I  — *  III,  I  — *■  IV)  are  shown  above 
the  curves. 

Results  and  Discussion 

Relative  Energies  of  Tetrazolium  Cations.  The  tetrazolium 
cation  has  four  parent  isomers  (all  R  =  H)  labeled  I,  II,  III,  and 
IV  (Figure  3).  Each  isomer  has  two  possible  resonance  contribu¬ 
tors.  The  relative  energies  of  I,  II,  III,  and  IV  are  shown  in 
Table  2.  Isomer  I  is  predicted  by  MP2  to  be  the  lowest  in  energy, 


with  isomer  II  only  1.9  kcal/mol  higher  in  energy.  The  MP2 
relative  energy  ordering  is  I  <  II  <  IV  <  III.  In  general,  DFT 
and  MP2  are  in  good  agreement,  although  DFT  slightly  reverses 
the  order  of  isomers  I  and  II. 

The  relative  energies  of  isomers  I  and  II  were  also  calculated 
using  the  CCSD(T)  method  at  the  MP2/6-31  lG(d,p)  geometries. 
Using  the  6-311G(d,p)  and  cc-pVTZ  basis  sets,  isomer  II  is 
predicted  to  be  lower  in  energy  than  I  by  0.2  and  0.5  kcal/mol, 
respectively.  The  open  chain  form  of  the  cation,  azidoformi¬ 
dinium  (V),  shown  in  Figure  4,  was  also  compared  to  the  four 
parent  isomers.  At  the  MP2/6-311++G(d,p)  geometry  this 
isomer  is  0.7  kcal/mol  lower  in  energy  than  isomer  I  (Table  3). 
According  to  CCSD(T)/cc-pVTZ//MP2/6-31+G(d),  isomer  V 
is  1.8  kcal/mol  higher  in  energy  than  isomer  I,  and  1.3  kcal/ 
mol  higher  in  energy  than  isomer  II. 

The  large  MP2/6-311G(d,p)  barriers  for  proton  transfer 
reactions  between  isomer  I  and  isomers  II,  III,  and  IV  (Figure 

5)  imply  that  movement  of  a  proton  on  the  ring  is  not  likely. 

Calculated  MP2/6-311G(d,p)  o  and  Jt  bond  orders40  (Figure 

6)  suggest  that  the  two  double  bonds  in  the  ring  are  delocalized. 
The  bond  lengths  between  atoms  in  the  rings  do  not  exhibit 
any  significant  changes  among  the  four  parent  isomers  (Table 
4).  The  bond  lengths  between  atoms  in  the  cation  rings  are 
consistent  with  a  delocalized  ring  system  (cf.,  1.32  A  for  a  CN 
double  bond  and  1.34  A  for  a  NN  double  bond). 

Bond  lengths  can  suggest  which  resonance  structures  are 
favored  for  isomers  I  and  III.  (The  two  resonance  structures 
for  isomers  II  and  IV  are  equivalent.)  The  nearly  equal  Nl  — 
N2  and  N2— N3  bond  lengths  in  isomer  I  (1.31  and  1.29  A, 
respectively)  suggest  that  resonance  structures  I  and  T  make 
similar  contributions  to  the  electronic  structure  of  the  ring,  with 
perhaps  a  slight  preference  for  I  over  T.  The  two  resonance 
structures  for  isomer  III  are  likewise  equally  important  as 
indicated  by  the  equal  N2— N3  and  N3— N4  bond  lengths  and 
the  nearly  equal  lengths  of  the  Nl— C5  and  N4— C5  bonds. 

MCSCF  jt  LMO  populations  and  Jt  bond  orders  (Table  5) 
can  provide  a  more  sophisticated  analysis  of  competing  reso¬ 
nance  structures.2  The  diagonal  density  matrix  elements  give 
the  electron  occupancy  of  the  localized  molecular  orbitals.  The 
off-diagonal  elements  give  the  bond  orders.  A  positive  bond 
order  indicates  a  bonding  interaction  and  a  negative  bond  order 
indicates  an  antibonding  interaction  between  the  LMOs. 


TABLE  2:  Relative  Energies  (kcal/mol)  of  Isomers  I,  II,  III,  and  IV 


cation 

B3LYP/6-31G(d,p) 

B3LYP/6-31  lG(d,p) 

MP2/6-31G(d,p) 

MP2/6-311G(d,p) 

I 

0 

0 

0 

0 

II 

-0.6 

-0.4 

1.4 

1.9 

III 

16 

16.2 

19.1 

19.2 

IV 

14.9 

14.7 

15.2 

14.7 

TABLE  3:  Relative  Energies  (kcal/mol)  of  the  Closed  Chain  Isomers  I  and  II  and  the  Open  Chain  Isomer  V,  Using  CCSD(T) 

MP2/6-3 1 1 ++G(d,p)  CCSD(T)/cc-pVTZ//MP2/6-3 1 1 ++G(d,p)  CCSD(T)/6-3 1 1  G(d,p)//MP2/6-3 1 1 ++G(d,p) 


I 

II 
V 


0.0 

1.9 

-0.7 


0.0 

-0.5 

1.8 


0.0 

-0.2 

1.6 


Tetrazolium-Based  Ionic  Liquids 


J.  Phys.  Chem.  B,  Vol.  110,  No.  23,  2006  11113 


a.) 


I-A 

II-A 

III-A 

IV- A 

H 

I 

H 

I 

H 

I 

H — C  N 

W  J/ 

N — N 

\ 

R 

H — C  N 

W+  // 

N — N 

/ 

R 

H — C  +N - R 

w  // 

N — N 

N 

H — QT  +N - H 

\  / 

N=N 

\ 

R 

b.) 

I-B 

II-B 

III-B 

IV-B 

H 

I 

H 

I 

H 

I 

R — cf  N 

W  J! 

N — N 

\ 

H 

R — C\  N 

W+  // 

N — N 

/ 

H 

R — C  +N - H 

w  // 

N — N 

N 

R — QT  +N - H 

\  / 

N=N 

\ 

H 

c.) 

I-C 

II-C 

III-C 

IV-C 

R 

I 

R 

I 

R 

I 

H — C  N 

W  J! 

N — N 

\ 

H 

H — C  N 

W+  // 

N — N 

/ 

H 

H — C  +N - H 

\\  // 

N — N 

.N 

H — +  N - R 

\  / 

N=N 

\ 

H 

Figure  7.  Isomers  with  a  single  substitution  on  a  nitrogen  or  a  carbon. 


TABLE  4:  Bond  Lengths  (A)  in  the  Ring  for  Isomers  I 
through  IV 


C5-N1 

N1-N2 

N2-N3 

N3-N4 

N4-C5 

I 

1.36 

1.31 

1.29 

1.33 

1.31 

II 

1.22 

1.35 

1.30 

1.35 

1.33 

III 

1.35 

1.32 

1.32 

1.32 

1.34 

IV 

1.35 

1.31 

1.32 

1.31 

1.35 

TABLE  5:  MCSCF 

it  Orbital  Populations  and  Bond  Orders 

71  orbital  populations 

Nl 

N2 

N3 

N4 

C5  total 

bonding 

I 

1.42 

1.10 

1.37 

1.12 

0.98 

1.93 

II 

1.44 

1.09 

1.09 

1.44 

0.94 

1.92 

III 

1.46 

1.46 

1.04 

1.10 

0.95 

1.91 

IV 

1.07 

1.43 

1.43 

1.07 

0.99 

1.94 

adjacent  jc  bond  orders 

C5-N1 

N1-N2 

N2-N3 

N3-N4 

N4-N5 

I 

0.55 

0.57 

0.65 

0.48 

0.73 

II 

0.65 

0.44 

0.77 

0.44 

0.65 

III 

0.61 

0.44 

0.60 

0.62 

0.64 

IV 

0.64 

0.56 

0.50 

0.56 

0.64 

next 

neighbor  jt  bond  ordrs  (antibonding) 

C5-N2 

N1-N3 

N2-N4 

N3-C5 

N4-N1 

I 

-0.22 

-0.30 

-0.27 

-0.11 

-0.15 

II 

-0.18 

-0.19 

-0.19 

-0.18 

-0.29 

III 

-0.30 

-0.33 

-0.22 

-0.15 

IV 

-0.14 

-0.34 

-0.34 

-0.14 

For  isomer  I,  the  Jt  orbital  populations  on  N1  and  N3  are 
1.42  and  1.37,  respectively.  In  a  resonance  structure  that  was 
purely  I,  N1  would  have  a  Jt  population  of  2.00  and  all  other 
atoms  in  the  ring  would  have  jt  populations  of  1 .00  (see  Figure 
3).  If  the  structure  were  purely  T,  then  N3  would  have  a  jt 
population  of  2.00  and  all  other  atoms  in  the  ring  would  have 


jt  populations  of  1.00  (cf.,  Figure  3).  In  reality  both  N1  and  N3 
have  Jt  populations  of  less  than  1.5  electrons.  This  indicates 
that  the  ring  is  a  hybrid  between  I  and  T  and  is  delocalized. 
The  difference  in  Jt  orbital  populations  indicate  that  I  is  slightly 
favored  over  T.  In  both  I  and  T  C5— N4  is  a  double  bond.  This 
bond  has  a  Jt  bond  order  of  0.73,  the  largest  of  the  five  bonds. 
The  N2— N3  Jt  bond  order,  a  double  bond  in  I,  is  0.65.  The 
Nl—  N2  jt  bond  order,  a  double  bond  in  T,  is  0.57.  Neither  of 
these  bond  orders  is  close  to  1.00  indicating  once  again  that 
the  ring  is  delocalized. 

A  similar  analysis  can  be  done  for  isomers  II— IV  showing 
that  all  the  Jt  electrons  are  delocalized  on  the  cations.  According 
to  the  Jt  orbital  populations  in  isomer  III,  the  lone  pair  is 
distributed  equally  on  Nl  and  N2.  Similar  comments  apply  to 
isomer  II(IV),  in  which  the  lone  pair  is  distributed  equally  on 
Nl  and  N4  (N2  and  N3). 

The  relative  energies  of  the  parent  isomers  with  a  single 
substituent  on  a  nitrogen  or  a  carbon  are  shown  in  Table  6. 
The  isomers  with  a  substituent  R'  (on  a  nitrogen)  are  labeled 
I- A  though  IV- A  (Figure  7a),  the  isomers  with  a  substituent 
R"  (on  the  carbon)  are  labeled  I-B  though  IV-B  (Figure  7b), 
and  the  isomers  with  a  substituent  R//r  (on  a  nitrogen)  are  labeled 
I-C  through  IV-C  (Figure  7c).  For  all  substituents,  the  relative 
energies  remain  the  same  as  those  of  the  parent  isomers  with 
no  substitution. 

A  — NO2  substituent  does  not  bind  well  to  a  nitrogen  of  the 
tetrazolium  ring:  the  NN  distance  is  2.4  A,  compared  with  a 
normal  single  NN  bond  distance  of  1.430  A  as,  for  example,  in 
N2H4.41  Because  of  the  lack  of  a  strong  NN  bond,  the  nitro 
group  easily  changes  positions  on  the  ring.  For  example,  III- 
A(NC>2)  converts  to  isomer  I-A(N02),  with  no  barrier.  A  cation 
with  a  — NO2  substitutent  bound  to  a  ring  nitrogen,  is  most  likely 
an  ion— dipole  complex,  in  which  the  cation  becomes  a  neutral 
tetrazole  ring  and  the  nitro  group  becomes  an  incipient  cation. 
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TABLE  6:  MP2/6-311G(d,p)  Relative  Energies  (kcal/mol) 
Isomers  with  a  Single  Substitution  on  a  Nitrogen  or  a 
Carbon" 


R' 

I— A 

II-A 

III-A 

IV- A 

H 

0.0 

1.9 

19.1 

14.7 

F 

0.0 

3.4 

18.3 

12.6 

CN 

0.0 

2.6 

18.6 

14.4 

n3 

0.0 

5.3 

16.8 

16.0 

nh2 

0.0 

1.8 

17.1 

13.0 

no2 

0.0 

2.7 

4.9 

R" 

I-B 

II— B 

III— B 

IV— B 

H 

0.0 

1.9 

19.1 

14.7 

F 

0.0 

5.0 

22.6 

16.8 

CN 

0.0 

1.3 

19.2 

14.6 

n3 

0.0 

2.9 

15.4 

11.0 

nh2 

0.0 

3.5 

16.3 

11.7 

no2 

0.0 

2.5 

17.8 

15.6 

R'" 

I-C 

II-C 

III-C 

IV-C 

H 

0.0 

1.9 

19.1 

14.7 

F 

0.0 

2.3 

22.1 

14.0 

CN 

0.0 

1.5 

19.9 

14.9 

n3 

0.0 

1.6 

20.1 

9.6 

nh2 

0.0 

1.0 

18.2 

9.4 

no2 

0.0 

1.1 

13.5 

3.3 

"  See  Figure  7  for  notation. 


TABLE  7:  MP2/6-311G(d,p)  Average  Energy  Differences 
(kcal/mol)  between  Substitution  at  N  and  C  on  the 
Tetrazolium  Ring 


average  values  (EC—  EN) 

F 

-56.9 

CN 

-20.1 

n3 

-30.0 

nh2 

-27.9 

no2 

+6.0 

This  is  easily  verified  by  examining  the  charge  on  the  substituent 
— NO2  group.  The  average  of  the  Mulliken  charges42  on  the 
— NO2  group  for  all  isomers  with  —  NO2  on  N  is  +0.83.  The 
average  MP2/6-31  lG(d,p)  binding  energy  of  NC>2+  to  a  tetrazole 
ring  is  22.8  kcal/mol. 

A  comparison  of  relative  energies  (averaged  for  all  four  parent 
isomers)  for  substitution  on  a  C  vs  N  is  shown  in  Table  7. 
Substitution  at  C  is  almost  always  energetically  favored  over 
substitution  at  N.  The  only  exception  is  for  the  —  NO2  substitu¬ 
ent.  The  bond  lengths  within  the  ring  do  not  change  significantly 
from  the  parent  isomer  when  a  single  substituent  is  present. 

Heats  of  Formation.  The  energy  content  of  a  heterocycle 
can  be  raised  dramatically  by  increasing  the  number  of  nitrogen 
atoms  in  the  ring.  For  example,  the  experimental  heats  of 
formation  for  imidazole,  1,2,4-triazole,  and  tetrazole  are  14.0, 
26.1,  and  56.7  kcal/mol,  respectively.49 

In  this  work,  heats  of  formation  for  the  cation  rings  were 
calculated  using  isodesmic  reactions  and  the  Gaussian-2  (G2) 
method,43’44  a  multilevel  method  designed  to  obtain  accurate 
thermochemistry.  The  G2  calculations  were  performed  using 
Gaussian  94  45  Isodesmic  reactions,  in  which  the  number  of 
formal  bond  types  is  conserved,  minimize  the  change  in 
correlation  energy,46  thereby  reducing  the  error  in  a  computed 
heat  of  formation. 

The  isodesmic  reactions  of  each  resonance  structure  are 
shown  in  Figure  8.  For  I  and  III  each  resonance  structure  has  a 
different  isodesmic  reaction.  The  heats  of  formation  for  I  and 
I'  (III  and  II T)  are  similar,  and  their  average  will  be  taken  as  a 
heat  of  formation  for  the  composite  I/T  (III/IIT).2  The  largest 


TABLE  8:  Heats  of  Formation  (kcal/mol),  from  G2  Theory, 
for  Parent  Isomers 


calcd  G3  heats 
of  formation" 

calcd  G2  heats 
of  formation^7 

I 

248.0 

245.1 

II 

247.1 

246.0 

III 

263.7 

264.3 

IV 

262.8 

259.8 

"  G3  Heats  of  Formation  are  taken  from  Satchell  and  Smith.19  b  From 
this  work. 


TABLE  9:  Changes  in  Heats  of  Formation  (kcal/mol),  from 
G2  Theory,  Due  to  Substituents,  Relative  to  Parent 
Compounds" 


R' 

I 

r 

av 

R" 

i  =  r 

R'" 

I 

I' 

av 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

17.1 

17.8 

17.5 

F 

-37.5 

F 

19.9 

20.6 

20.2 

CN 

70.7 

71.7 

71.2 

CN 

51.0 

CN 

70.5 

71.5 

71.0 

n3 

108.6 

111.4 

110.0 

n3 

79.7 

n3 

109.9 

112.8 

111.3 

nh2 

18.2 

18.9 

18.6 

nh2 

-10.0 

nh2 

20.7 

21.4 

21.0 

no2 

31.1 

18.1 

24.6 

no2 

18.0 

no2 

32.8 

19.7 

26.3 

R' 

II 

ir 

av 

R" 

ii  =  ir 

R'" 

II 

II' 

av 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

20.3 

21.0 

20.6 

F 

-36.0 

F 

20.3 

21.0 

20.7 

CN 

70.2 

71.2 

70.7 

CN 

51.7 

CN 

70.1 

71.1 

70.6 

n3 

109.6 

112.5 

111.0 

n3 

83.1 

n3 

109.6 

112.5 

111.1 

nh2 

19.9 

20.6 

20.3 

nh2 

-10.0 

nh2 

19.8 

20.5 

20.1 

no2 

32.0 

18.9 

25.4 

no2 

18.6 

no2 

32.0 

18.9 

25.4 

R' 

III 

nr 

av 

R" 

hi  =  nr 

R'" 

III 

III' 

av 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

20.5 

21.2 

20.8 

F 

-38.4 

F 

22.8 

23.5 

23.1 

CN 

70.8 

71.8 

71.3 

CN 

50.3 

CN 

71.2 

72.3 

71.8 

n3 

104.8 

107.6 

106.2 

n3 

77.3 

n3 

110.8 

113.7 

112.2 

nh2 

15.4 

16.0 

15.7 

nh2 

-12.1 

nh2 

19.8 

20.4 

20.1 

no2 

no2 

16.6 

no2 

27.1 

14.0 

20.6 

R' 

IV 

IV' 

av 

R" 

IV  =  IV' 

R'" 

IV 

IV' 

av 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

19.2 

19.9 

19.5 

F 

-39.6 

F 

19.2 

19.9 

19.6 

CN 

70.7 

71.7 

71.2 

CN 

50.7 

CN 

70.7 

71.7 

71.2 

n3 

104.9 

107.8 

106.3 

n3 

81.0 

n3 

104.8 

107.7 

106.3 

nh2 

15.3 

15.9 

15.6 

nh2 

-11.7 

nh2 

15.4 

16.1 

15.7 

no2 

21.3 

8.3 

14.8 

no2 

18.9 

no2 

21.4 

8.3 

14.9 

Refer  to  Figure  3  for  structures. 


difference  in  heats  of  formation  between  resonance  structures 
is  13.0  kcal/mol  for  the  —  NO2  substituent;  the  other  differences 
are  all  less  than  3  kcal/mol. 

The  heats  of  formation  for  the  parent  isomers  are  shown  in 
column  2  of  Table  8.  These  compare  well  to  heats  of  formation 
that  were  calculated  previously  using  Gaussian-3.19  The  changes 
in  heats  of  formation  for  substituted  cations  relative  to  the 
hydrogen  substituted  cations  are  shown  in  Table  9.  For  the  most 
part,  substitution  destabilizes  the  ring,  although,  substitution  by 
a  fluorine  or  an  amine  on  a  carbon  stabilizes  the  ring  by  up  to 
40  kcal/mol.  However,  an  azide  group  or  a  nitrile  group  can 
increase  the  heat  of  formation  of  the  ring  by  as  much  as  112 
kcal/mol.  Substitution  with  — F,  —  NH2,  and  —  NO2  generally 
increases  the  heats  of  formation  by  about  20  kcal/mol. 

Proton-Transfer  Reactions.  One  possible  reaction  pathway 
for  the  simple  ion  pairs  is  proton  transfer  from  the  cation  to  a 
partner  anion  to  form  a  neutral  pair.  The  heats  for  the  reaction 
CN4H3+  +  y-  —  CN4H2  +  HY  (where  Y  =  N03,  C104,  or 
N(N02)2)  are  shown  in  Table  10.  All  of  the  energy  differences 
(A E)  are  exothermic,  by  46—144  kcal/mol,  and  the  AEs  for 
proton  transfer  to  each  of  the  three  anions  are  similar.  The  heat 
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Figure  8.  Isodesmic  reactions  for  all  resonance  structures. 


TABLE  10:  Heats  of  Proton  Transfer  Reaction  (kcal/mol) 
for  Anion  =  N03_,  C104~,  and  N(N02)2~  for  Deprotonation 
at  H',  H",  and  HI"  (See  Figure  3) 


cation 

anion 

H' 

H" 

Hi" 

I 

NO3- 

-124.8 

-86.1 

-128.5 

cio4- 

-102.6 

-63.9 

-106.3 

N(N02)2- 

-112.0 

-73.3 

-115.7 

II 

no3- 

-126.6 

-103.7 

same  as  H' 

CKV 

-104.4 

-81.5 

same  as  H' 

N(N02)2- 

-113.9 

-90.9 

same  as  H' 

III 

no3- 

-143.7 

-85.9 

-147.4 

CKV 

-121.5 

-63.7 

-125.3 

N(N02)2- 

-131.0 

-73.1 

-134.7 

IV 

no3- 

same  as  Hi" 

-68.9 

-142.9 

C1(V 

same  as  H"' 

-46.7 

-120.8 

N(N02)2- 

same  as  Hi" 

-56.2 

-130.2 

of  reaction  for  removal  of  Hi" ,  shown  in  column  three  of  Table 
10,  is  more  exothermic  than  the  removal  of  H'  (see  Figure  3) 
because  the  2,5-disubstituted  neutral  ring  is  3.7  kcal/mol  lower 
in  energy  than  the  1,5-disubstitued  neutral  ring.  The  1,3- 
disubstituted  ring  is  considerably  higher  in  energy.  There  are 


H2N-R’  +  H3N+-R’”  +  H3C-R” 
+  2  H2N-NH2  +  H2C=NH  + 
HN=NH  +H3C-NH2 


three  possible  positions  for  proton  transfer  to  the  dinitramide 
anion,  but  dinitramic  acid  is  most  stable  with  the  proton  on  the 
center  nitrogen,  so  this  is  the  acid  used  when  calculating  the 
AEs. 

Ion  Pair  Interactions.  Ionic  liquid  dimer  pairs  were  opti¬ 
mized  with  DFT  and  MP2  methods  using  the  6-31+G(d)  basis 
set.  ZPE  corrections  at  the  B3LYP/6-31+G(d)  level  of  theory 
were  used  in  all  cases.  The  structures  of  the  tetrazolium  dimers 
with  perchlorate  and  nitrate  partners  are  shown  in  Figures  9 
and  10,  respectively.  For  the  cation  with  perchlorate  or  nitrate, 
at  both  levels  of  theory  a  proton  transfers  from  the  cation  to 
the  anion  during  geometry  optimization  of  the  ion  pair. 
Hydrogen  bonding  and  other  attractive  electrostatic  interactions 
are  found  to  be  important  in  the  dimer  structures.  Optimizations 
give  geometries  in  which  the  number  of  attractive  electrostatic 
interactions  is  maximized.  For  pairs  containing  HNO3  or  HCIO4, 
there  is  always  one  linear  O— H — N  hydrogen  bond.  Hydrogen 
bond  lengths  vary  from  1.795  to  1.846  A  for  HNO3  and  1.771 
A  to  1.822  A  for  HCIO4.  The  majority  of  structures  also  contain 

a  nonlinear  N— H - O  hydrogen  bond  or  a  C— H - O  attractive 

electrostatic  interaction. 
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Figure  9.  Optimized  structures  (A)  and  MP2  and  B3LYP  (in 
parentheses)  relative  energies  (kcal/mol)  for  proton  transferred  structures 
that  result  from  optimization  of  initial  geometries  of  a  tetrazolium  cation 
paired  with  a  NC>3~  anion.  MP2  heats  of  reaction  (in  brackets)  in  kcal / 
mol  for  the  following  reaction:  tetrazolium  +  NO3-  — *  [tetrazole]  [H— 
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Figure  10.  Optimized  structures  (A)  and  MP2  and  B3LYP  (in 
parentheses)  relative  energies  (kcal/mol)  for  proton  transferred  structures 
that  result  from  optimization  of  initial  geometries  of  a  tetrazolium  cation 
paired  with  a  CIO4-  anion.  MP2  heats  of  reaction  (in  brackets)  in  kcal/ 
mol  for  the  following  reaction:  tetrazolium  +  CHAT  [tetrazole]  - 
[H-CIO4]. 


Tetrazolium  dinitramide  (Figure  11)  pairs  have  a  much  wider 
variety  of  structures  than  tetrazolium  nitrate  or  tetrazolium 
perchlorate  pairs.  This  is  because  the  dinitramide  anion  has  the 


ability  to  bind  a  proton  at  three  different  positions  (Figure  12). 
For  HN(N02)2  pairs,  hydrogen  bond  lengths  vary  from  1.721 
to  1 .942  A.  Once  again  hydrogen  bonding  and  other  electrostatic 
interactions  are  important  for  structures  containing  HN(N02)2. 
However,  6  out  of  19  structures  do  not  contain  a  linear  hydrogen 
bond  at  all.  These  structures  usually  contain  weaker,  nonlinear, 

N— H - N,  O— H - N,  or  N— H - O  attractive  electrostatic 

interactions. 

The  relative  energies  of  the  pairs  with  HNO3  or  HCIO4  show 
that  the  lowest  energy  pairs  contain  the  1,3-H  substituted 
tetrazole  isomer.  For  example,  structure  1  in  Figure  11  is  2.8 
kcal/mol  lower  in  energy  than  structure  6.  This  difference  in 
energy  can  be  attributed  to  the  difference  in  energy  of  the  1,3-H 
substituted  tetrazole  ring  (number  1)  and  the  1,2-H  substituted 
tetrazole  ring  (number  6).  The  lowest  energy  structures  with 
HN(N02)2,  contain  the  1,3-H  substituted  tetrazole  isomer  and 
have  the  proton  transferred  to  the  central  nitrogen  of  N(N02)2_. 
Also  shown  in  Figures  9—11  are  the  changes  in  energy  for  the 
reaction  in  which  the  infinitely  separated  anion  and  cation  come 
together  to  form  either  an  ion  pair  or  a  neutral  pair  by  proton 
transfer  from  the  cation  to  the  anion: 

[tetrazolium] +  +  [anion]-  — *  [tetrazole] [H- anion] 


or 

[tetrazolium] +  +  [anion]  -  — *  [tetrazolium]  [anion] 
(anion  =  N03-,  C104-,  or  N(N02)2-) 

For  pairs  in  which  a  proton  transfers,  the  neutral  pair  is  lower 
in  energy  than  the  separated  neutral  products,  due  to  the 
stabilization  from  electrostatic  interactions  in  the  neutral  dimer 
products.  The  average  dimer  stabilization  is  7.3,  9.6,  and  8.4 
kcal/mol  for  [HN03],  [HCIO4],  and  [HN(N02)2]  respectively. 

At  the  DFT  level  of  theory  there  are  some  stable  N(N02)2~— 
cation  pair  geometries;  however  when  those  DFT  geometries 
are  optimized  with  MP2,  the  acidic  proton  usually  transfers  to 
the  anion.  Only  one  stable  anion— cation  pair  was  found  at  the 
MP2  level  of  theory.  This  ion  pair  is  structure  13  in  Figure  11. 
The  ion  pair  is  higher  in  energy  by  17.8  kcal/mol  than  the  lowest 
energy  neutral  pair  at  MP2/6-31+G(d).  In  this  structure  the 
N— H  bond  length  is  stretched  to  1.076  A  as  compared  to  1.019 
A  for  the  nonacidic  hydrogen.  The  hydrogen  bond  length  is 
also  much  shorter  at  1.598  A.  The  length  of  the  hydrogen  bonds 
show  that  the  hydrogen  bond  in  the  ionic  structure  is  much 
stronger  than  the  hydrogen  bond  in  the  various  neutral  stmctures. 
The  MP2  barrier  for  proton  transfer  from  the  cation  to  the  anion 
of  structure  13  of  Figure  11  is  0.4  kcal/mol  (Figure  13); 
however,  when  ZPE  corrections  are  accounted  for,  the  barrier 
disappears.  For  all  dimer  pairs,  proton  transfer  from  a  cation  to 
the  nitrogen  of  the  dinitramide  anion  is  favored  energetically 
over  proton  transfer  to  one  of  the  oxygens  of  the  dinitramide 
anion. 

Gas  phase  calculations  on  an  ionic  pair  can,  of  course,  only 
give  an  approximation  to  the  true  interactions  in  a  crystallized 
ionic  liquid.  Proton  transfer  from  the  cation  to  the  anion  should 
be  less  likely  in  the  crystal  or  in  bulk  liquid.  For  example,  in 
studies  of  ammonium  salts,  proton  transfer  occurs  in  the  isolated 
gas-phase  ion  pair,  but  stable  ion  pairs  are  found  when  two 
anions  and  two  cations  are  present.47  MP2/6-31++G(d,p) 
calculations  on  1,2,4-triazolium  dinitramide  typically  do  not  give 
stable  ion  pairs.2  The  most  recent  MP2  calculations  done  on 
three  ion  pairs  of  1,2,4-triazolium  dinitramide,  show  that  the 
six  ion  cluster  is  only  1.0  kcal/mol  higher  in  energy  than  the 
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Figure  11.  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative  energies  (kcal/mol)  of  the  tetrazolium  cation  paired  with  a 
N(N02)2-  anion.  MP2  heats  of  reaction  (in  brackets)  in  kcal/mol  for  the  following  reaction:  tetrazolium  +  N(N02)2-  [tetrazole][HN(N02)2]  or 
(see  text)  [tetrazolium]  [N(N02)2]. 
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Figure  12.  Three  different  geometries  of  dinitramic  acid. 


I* 


Figure  13.  MP2/6-31+G(d)  minimum  energy  path  for  proton  transfer 
from  the  1,4,5-tetrazolium  cation  to  the  dinitramide  anion  to  form  1,5- 
tetrazole  and  dinitramic  acid.  The  relative  energies  including  ZPE 
corrections  are  in  parentheses.  Relative  energies  are  in  kcal/mol. 

six  neutral  cluster.48  This  is  caused  by  an  increase  in  charge 
balance  that  occurs  as  the  number  of  ion  pairs  present  increases. 
That  is,  each  negative  (positive)  ion  is  balanced  by  an  increasing 
number  of  positive  (negative)  ions.  The  structures  and  relative 
energies  for  multiple  ion  pairs  is  still  under  investigation.  Future 
studies  on  larger  clusters  of  ionic  pairs  are  needed  to  more  fully 
understand  the  interactions  in  the  bulk  liquid. 

Conclusion 

In  this  work,  tetrazolium  cations  I  and  II  were  established  to 
be  the  lowest  energy  isomers,  for  the  parent  isomers  as  well  as 
for  all  substituted  isomers.  The  DFT  and  MP2  energies  are 
generally  in  good  agreement.  The  relative  energy  of  the  open 
chain  form  of  the  cation  is  predicted  to  be  only  slightly  higher 


in  energy  than  isomers  I  and  II.  A  MCSCF  Jt  orbital  analysis 
indicates  that  the  electrons  in  the  cation  ring  are  delocalized. 
Calculated  heats  of  formation  show  that  the  tetrazolium  cation 
ring  has  the  potential  to  release  large  amounts  of  energy  during 
decomposition  and  thus  has  excellent  potential  as  a  high  energy 
fuel.  This  is  especially  true  when  the  ring  is  substituted  with 
— N3  or  — CN.  When  a  cation  is  paired  with  oxygen  rich  anions, 
a  single  gas  phase  ion  pair  was  not  generally  found  to  be  stable. 
A  proton  transfers  without  barrier  from  the  cation  to  the  anion 
to  form  a  neutral  pair. 
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